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On the ergodicity bounds for a constant retrial rate queueing model 
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Abstract. We consider a Markovian single-server retrial queueing system with a constant retrial 
rate. Conditions of null ergodicity and exponential ergodicity for the corresponding process, as well as 
bounds on the rate of convergence are obtained. 

Key words: single-server retrial queueing system; constant retrial rate; continuous-time Markov 
chains; ergodicity bounds 

1. Introduction 

We consider the following Markovian single-server retrial queueing system with a constant 
retrial rate, denoted by S. The exogenous (primary) customers follow a Poisson input with rate 
A. The customers have i.i.d. exponential service times {S^}, with a generic element S and rate 
/i := l/ES*. If a new customer Ends a server busy it joins an inhnite-capacity orbit and called 
secondary customer. We assume that the orbit works as a single FIFS (£rst-in, hrst-served or 
first-come, hrst-served) server. It means that, if the orbit is not empty, a head line (the oldest) 
secondary customer attempts to enter the server after an exponentially distributed time with a 
rate po- Thus, unlike classical retrial models, the orbit rate in S does not depend on the orbit 
size, i.e., the number of secondary customers. Such a model is called a retrial model with constant 
retrial rate. Thus, orbit can be interpreted as a single-server -/M/l-type queue with service rate 
Po, and input is generated by the customers rejected in busy server. Note that the only possible 
source of instability of such system is an inhnite growth of the orbit size. A sufficient stability 
condition of the general single-class retrial system with constant retrial rate described above is 
obtained in [S]. Stability analysis of multi-class retrial system is presented in [B]. 

Considered system can be successfully applied to model the standard multi-access protocol 
ALOHA, with restrictions for the individual retrial rates and has other several applications. We 
list the most important papers related to motivation of presented model. In |T0] Fayolle first used 
a retrial queue with constant retrial rate to simulate a telephone exchange system. Then, in 
the authors have modelled an unslotted Carrier Sense Multiple Access with Collision Detection 
(CSMA/CD) protocol and in [8] and [9] the authors have modelled some particular versions of the 
ALOHA protocol. The model of Fayolle [10] was extended in n. PI for more complex settings, such 
as multiple servers and waiting places. In [T2| it has been proposed to use the retrial queue with 
constant retrial rate to model a logistic system. In [3] and |1| the authors have suggested to use 
retrial queues and retrial networks with constant retrial rates to model TCP traffic originated from 
short HTTP connections. The present retrial model also appears to be relevant for the optical- 
electrical hybrid contention resolution scheme for Optical Packet Switching networks Hal HD. 

In this note we hrstly consider the description of the model (Section 2). In Section 3, we present 
the auxiliary results that were obtained for birth-death processes in jlTl [151 HB] • In Sections 4, 
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5 we apply this approach for the considered model and obtain explicit bounds on the rate of 
convergence both for null ergodic and exponentially ergodic situations. 

2. Description of the model 

Now we describe the model in more detail. For a moment t, let z/(t) be the number of customers 
in the server and N{t) be the number of customers in the orbit. That is, u{t) = 0, if the server is 
empty, and z/(t) = 1, otherwise, while N{t) = 0, 1, .... We introduce the basic two-dimensional 
process X(t) = N(t), t > 0} with the state space {0, 1} x {0, 1, 2, ...}. Consider in more 

detail the transitions between the states of the system. To this end, we first enumerate the states 
of the process as follows: each state {0, n} will be denoted 2n -|- 1, for n > 0, while each state 
{1, n} is denoted 2n, n > 1. Denote hj Q = {qij) an intensity matrix corresponding to the given 
enumeration. Then it follows that 


and, for n > 1, 


gi,2 = A, 


Q2n,2n 

Q^2n,2n+2 

q2n,2n-l 

Q2n-\-l,2n-\-l 

Q^2n+l,2n 

Q'2n+l,2n+2 


—(A -I- /i), 

A, 

/X, 

— (A -I- /io), 
hO; 

A. 


As a result. 

an intensity 

matrix Q = (qtj) 

takes the form 






/-A 

A 

0 

0 

0 

0 

0 

0 

0 


fi 

-(A + 

h) 0 

A 

0 

0 

0 

0 

0 


0 

ho 

— (A -|- /io) 

A 

0 

0 

0 

0 

0 

Q = 

0 

0 

h 

—(A -f /i) 

0 

A 

0 

0 

0 

0 

0 

0 

ho 

—(A -|- /io) 

A 

0 

0 

0 


0 

0 

0 

0 

h 

— (A -f /i) 

0 

A 

0 


v ^ 










Figures 1,2 illustrate these transitions by different but equivalent ways: 
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3. Auxiliary notions 


Then the probabilistic 
system: 


dynamics of the process is represented by the forward Kolmogorov 


dp 
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Figure 1: 



where A = is the corresponding transposed intensity matrix, and 

p = p{t) = (piit),p 2 it), ... )^ is the column vector of state probabilities of the process X{t). 

Throughout the paper by || ■ || we denote the Zi-norm, i. e. ||x|| = ^ |xi|, and ||i?|| = 
suPj \bij\ for matrix B = (%)“•=!. 

We recall here the approach for the study on the rate of convergence of birth-death process, 
see all details and discussion in [TTl [151 (13 HE] • 

Let 5 be a bounded linear operator on a Banach space B and let / denote the identity operator. 
The number 
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is called the logarithmic norm of B. 

\i B = /i, so that the operator B is given by the matrix B = then the logarithmic 

norm of B can be found explicitly: 
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On the other hand, the logarithmic norm of the operator B is related to the Cauchy operator 


3 















V(t, s) of differential equation 


in the following way 
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t > 0. 


h —^-j-O h 

From the latter one can deduce the following bounds of the Cauchy operator V(t,s), see 

[IIIII51II6]: 


l|V'(i,s)llB 0<s<i, 

Let n be a set of all stochastic vectors, i. e. li vectors with nonnegative coordinates and unit 
norm. Hence ||H|| = 2A + 2max(/i,/ xq) < oo. Thus, the Cauchy problem for differential equation 
Q has a unique solutions for an arbitrary initial condition, and p(s) G implies p(t) G for 
f > s > 0. 


4. Null ergodicity 

We recall that a Markov chain X{t) is called null ergodic (in other terms, it means null 
recurrence or transience), if Pk{t) —)■ 0 as t —)■ oo for any initial condition p(0) and any k. 

Consider a decreasing sequence of positive numbers i = 1,2,..., and the corresponding 
diagonal matrix A with diagonal entries { 6 k}- 
Let ZiA be the space of sequences: 


hA = {p = {Pi,P2, • • • : IIpIIia = ||Ap|| < oo} . 


Then we have 





( 2 ) 


Put (5i = 1 and 62 k+i = 0 ^ 2 /^ 62 k = b 62 k-i, if /c > 1 for some a,b. Then 62 k+i = 
62 k = a^~^b^, if fc > 1. 
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A(1 — 6) + /xo ’ b\ 

Then a < 1,6 < 1, A(1 — ab) — ^{b~^ — 1) > 0, A(1 — b) — ho(o~^ — 1) > 0. 

Hence ineqnality < 0 implies nnll ergodicity of X{t), and we obtain the following 

statement. 


Theorem 1. Let assumption Q) hold. Then the process X{t) is null ergodic and 
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for any t > 0, any initial condition X(0) = k, and any natural N, where 

Q* = min (A(l — ab) — ia{b~^ — 1), A(1 — b) — ho(a~^ — 1)) > 0. 
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5. Exponential ergodicity 


By introducing pi (t) 


1 — from (jT| we obtain the equation 
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see detailed discussion in dg ng [ig. Let {dj}, i = 2, 3,..., be a sequence of positive numbers. 
Put Pi = Ylh =2 ^n- 

Let D be the upper triangular matrix. 
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and liD be the corresponding space of sequences 
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where 


0^2 — X fj^ — A (6 (ife), 
a 2 fe+i = A +/io -^>1, 

C( 2 k +2 = A + /i — A(fe + ab) — fioa~^, k > 1. 

Hence exponential ergodicity of the process X{t) follows from the bonnd 

inf tti = min (A + po — A + /i — X{b + ab) — fXoa~^) > 0, 


for some a, b snch that ab > 1 , see HH Ea ES]. 
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statement. 




Theorem 2. Let assumption (15) be true. Then the process X{t) is exponentially ergodic and the 
following bound holds: 

~ (19) 
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for any t > 0 , and any initial condition p( 0 ), where 

a* = min (A + po — A + p — X{b + ab) — kooT^) > 0, 
and 71 = (ttojTTi, ... )'^ is the corresponding stationary distribution. 


Remark 1. It is worth mentioning that condition (15) is the stability criteria of the system studied 
by regenerative method in 13/. More exactly, under this condition, the workload and gueue size 
processes are positive recurrent regenerative (that is, the mean regeneration period is finite). 


Remark 2. One can obtain the respective perturbation bounds applying Theorem 2 and approach 
ofinS' On the other hand, the opposite condition iB implies an unlimited growth of the process, 
and it is consistent with the null ergodicity which established in Theorem 1. 
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